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We study the magnetic properties of nanoscale magnetic 
films with large perpendicular anisotropy comparing polar- 
ization microscopy measurements on Co2sPt72 alloy samples 
based on the magneto-optical Kerr effect with Monte Carlo 
simulations of a corresponding micromagnetic model. In our 
model the magnetic film is described in terms of single-domain 
magnetic grains, interacting via exchange as well as via dipo- 
lar forces. Additionally, the model contains an energy barrier 
which has to be overcome in order to reverse a single cell and 
a coupling to an external magnetic field. Disorder is taken 
into account. 

We focus on the understanding of the dynamics especially 
the temperature and field dependence of the magnetisation 
reversal process. The experimental and simulational results 
for hysteresis, the reversal mechanism, domain configurations 
during the reversal, and the time dependence of the magneti- 
sation are in very good qualitative agreement. The results 
for the field and temperature dependence of the domain wall 
velocity suggest that for thin films the hysteresis can be de- 
scribed as a depinning transition of the domain walls rounded 
by thermal activation for finite temperatures. 
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I. INTRODUCTION 

In recent years great effort was focused on the mag- 
netisation reversal process in magnetic thin films with 
perpendicular anisotropy because of their potential ap- 
plication as high density recording media. In particular 
CoPt alloy films were found to be a very promising com- 
pound for magnetic and magneto-optic storage jlj. 

Two different mechanisms can be thought of to dom- 
inate the reversal process: either nucleation or domain 
wall motion Which of these mechanisms dominates 
a reversal process depends on the interplay of the dif- 
ferent interaction forces between domains with different 
magnetic orientation. In recent experiments on Co2sPt72 
alloy films |3|Q] a crossover from magnetisation reversal 
dominated by domain growth to a reversal dominated by 
a continuous nucleation of domains was found depend- 
ing on the film thickness which was varied from lOnm to 



30nm. Correspondingly, characteristic differences for the 
hysteresis loops have been found. Similar results have 
been achieved by simulations of a micromagnetic model 
using zero temperature dynamics |3|Q| and Monte Carlo 
methods respectively 0. 

It is the goal of this paper to work out the relation 
between experiments and simulations. An exact, quanti- 
tative description of the experimental results is hindered 
by the facts that the observable length scales in exper- 
iment and simulation are different and that there is no 
straight forward mapping of the time scale of a Monte 
Carlo simulation on experimental time scales. However, 
the emphasis of our work is on a deeper qualitative under- 
standing of the dynamical aspects of the magnetization 
reversal and especially on the principal influence of the 
field and the temperature on the domain wall velocity. 
We find that for thin films these dynamics is governed 
by a depinning transition of the domain walls rounded 
by thermal activation for finite temperatures. 



II. EXPERIMENTAL METHODS 

The CoPt alloy films were prepared under UHV condi- 
tions at low deposition rates on Si substrates covered with 
a 20nm Pt buffer layer. In order to achieve a dominant 
perpendicular anisotropy the films were grown at rather 
high substrate temperatures of about 220°C. The compo- 
sition of 28 at.% Co and 72 at.% Pt is known to exhibit 
the maximum of the polar Kerr rotation for the whole 
range of alloys pL Structural characterization reveals 
a predominant fcc(lll) texture. The films were usually 
found to form a polycrystalline microstructure. Scanning 
tunneling microscopy was used to determine the average 
grain size. The preparation process yields rather uniform 
grain size of the order of 20nm almost independent of the 
film thicknesses, which cover the range of 5 to 50nm. 

X-ray analysis reveals a dispersion of the crystalline 
c-axis in the range of about 3 degrees. For the mag- 
netic characterization we applied magneto-optic Kerr mi- 
croscopy at very high spatial resolution of about 1 fim [Q . 
The typical Kerr rotation for the CoPt alloy films is rang- 
ing from 0.1° to 0.3°. High speed image processing tech- 
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nique using a low noise CCD camera is able to acquire 
domain patterns with millisecond time resolution. The 
experimental setup permits to gain complete hysteresis 
loops with applied external magnetic field up to 0.5T. 
An electrical heater assembly provides the variation of 
the sample temperature by direct heating using a bifi- 
lar heat wire. Image processing software was developed 
to determine the domain wall motion and the average 
magnetisation of the sample. 

III. MICROMAGNETIC MODEL 

Co28Pt72 alloy films have a polycrystalline structure. 
For a theoretical description by a micromagnetic model 
H the film is thought to consist of cells on a square lattice 
with a square base of size L x L where L = 20nm. The 
height h of the cells is varied from lOnm to 30nm. Due 
to the large anisotropy of the CoPt alloy film the cells 
are thought to be magnetised perpendicular to the film 
only with a uniform magnetisation M s which is set to the 
experimental value of M s = 365kA/m for the saturation 
magnetisation in these systems . The cells interact via 
domain wall energy and dipole interaction. The coupling 
of the magnetisation to an external magnetic field H is 
taken into account as well as an energy barrier which has 
to be overcome during the reversal process of a single 
cell. 

From these considerations it follows that the change of 
energy caused by reversal of a cell i with magnetisation 
L 2 hM s Oi with Oi = ±1 and Act = a(new) — a(old) = ±2 
is: 

AEi = -^LhS w A<7i a j 

(3) 

+^M!Lh 2 Aa i Y^v(a j ,r ij ) 
3 

—[i HL 2 hM s A(Ji (1) 

The first term describes the wall energy AE W and the 
sum is over the four next neighbors. One can expect 
that the grain interaction energy density S w is a reduced 
Bloch-wall energy density Sb since the crystalline struc- 
ture of the system is interrupted at the grain boundary 
and also since there may be a chemical modulation of the 
CoPt alloy at the grain boundary. We varied the value of 
S w in the simulation and got the best agreement with ex- 
perimental results for a value of S w — 0.0022J/m 2 which 
is approximately 50% of the Bloch-wall energy Sb for 
this system 0. 

In the second term describing the dipole coupling AEd 
the sum is over all cells, ry is the distance between two 
cells i and j in units of the lattice constant L. For large 
distances it is v(<jj,rij) = -^f . For shorter distances a 

more complicated form which is a better approximation 
for the shape of the cells which we consider can be deter- 
mined numerically and was taken into account. 



The third term describes the coupling AEh to an ex- 
ternal field H. 

Additionally, an energy barrier Si must be considered 
which describes the fact that a certain energy is needed 
to reverse an isolated cell. Two reversal mechanisms can 
be taken into account as limiting cases: (i) coherent ro- 
tation of the magnetisation vector described by an angel 
9: In this case the anisotropy leads to an energy barrier 
of L 2 hK u where K u is the anisotropy constant which is 
K u — 200k J/m 3 0. (ii) domain wall motion through 
the grain: In this case the energy barrier is LhSb due 
to the fact that the domain wall energy is lowered at 
the grain boundary. The highest possible value for Sb 
is the Bloch-wall energy mentioned above but it can be 
assumed that the for the energy barrier relevant value 
is smaller than the value above since a Bloch wall is 
thicker than the size of the cell. Comparing these two 
energies for the reversal mechanism of a CoPt grain one 
finds that here domain wall motion through the grains 
has the lower energy barrier so that from now on only 
this mechanism will be taken into account. We assume 
that during the reversal process the energy barrier has 
its maximum value LhSb when the domain wall is in the 
center of the cell, i.e. when half of the cell is already 
reversed. Consequently, the energy barrier which has to 
be considered in a Monte Carlo simulation is reduced to 
S = max(0, LhS b ~ \\{E W + E d + E H )\). The simula- 
tions are in good agreement with the experiments using 
S b = 0.0007J/m 2 . 

Note, that we have introduced two different wall energy 
densities, S w for the nearest neighbor interaction, and 
Sb for the intrinsic energy barrier. They have the same 
physical origin but they have to be handled separately in 
our model: e. g., the limit of a system of isolated grains 
is described by S w — ► but then there is still a finite Sb 
due to the existence of energy barriers which are relevant 
for the flip of the isolated cells. The opposite limit is a 
system without grain boundaries. Here it is Sb = 0, i. e., 
there is no energy barrier but the domain wall has still 
a constant, finite energy density S w which in this limit 
does not depend on the position of the wall. 

Obviously, the grain sizes in the CoPt alloy are ran- 
domly distributed ||. In order to be realistic, in our 
simulation disorder has to be considered (see also |8| for 
a discussion of the role of disorder in a similar simula- 
tion). In the model above this corresponds to a random 
distribution of L which can hardly be simulated exactly 
since it modulates the normalized cell distance of the 
dipole interaction. Therefore, as a simplified ansatz to 
simulate the influence of disorder we randomly distribute 
L in the energy term that describes the coupling to the 
external field. Here a random fluctuation of L is most 
relevant, since this term is the only one scaling quadrat- 
ically with L. In the simulations we use a distribution 
which is Gaussian with width A = 0.1. Through this 
kind of disorder our model is mapped on a random-field 
model. Other possible origins of disorder in the CoPt 
films are e. g. fluctuations of the easy axis of the grains 
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or fluctuations of the saturation magnetisation. Interest- 
ingly these kinds of disorder would also predominantly 
influence the last term in Eq. |l| and hence can be gath- 
ered in the random field. 

The simulation of the model above was done as in ear- 
lier publications |||J via Monte Carlo methods || using 
the Metropolis algorithm with an additional energy bar- 
rier. Since the algorithm satisfies detailed balance and 
Glauber dynamics it allows the investigation of thermal 
properties as well as the investigation of the dynamics of 
the system. Note that in the limit of low temperatures 
the Monte Carlo algorithm passes into a simple energy 
minimization algorithm with single spin flip dynamics, 
so that also the case of zero temperature can be investi- 
gated. 

The size of the lattice was typically 150 x 150. The 
dipole interaction was taken into account without any 
cut-off or mean field approximation. 

IV. HYSTERESIS AND THE REVERSAL 
MECHANISM 

We start our analysis by comparing the hysteresis loop 
from simulations for T — 300K, Fig. with the corre- 
sponding experimental hysteresis loop for room temper- 
ature, see Ref. [|]|| . 
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FIG. 1. Simulated hysteresis loops for a 30nm and a lOnm 
film at T = 300K. 



The qualitative agreement is very good. In both cases, 
experiment as well as simulation the hysteresis loop of 
the thin film is nearly rectangular. For the thicker film 
the loop has a finite slope, especially at the end of the 
hysteresis. This can be understood through the different 
reversal mechanisms. The reversal of the lOnm film is 
dominated by domain wall motion. Once a nucleus be- 
gins to grow the domain wall motion does not stop until 
the magnetisation has completely changed. For the 30nm 
film the enhanced dipolar forces stabilize a mixed phase 
which can be changed only by a further increase of the 
external field §-§. The reason for the enhanced influ- 
ence of the dipolar forces can be seen in Eq.|l| where the 
dipole term is the only term that scales quadraticly with 
the film thickness h which means that the influence of the 
dipole coupling is neglectable for h — > and on the other 
hand strongly growing with increasing film thickness. 

Quantitatively, the agreement of the simulation with 
the experimental results is reasonable only for the 30nm 
film while the nucleation field is much to high for the 
simulation of the lOnm film. There are several possible 
reasons for this effect. First, the nucleation field which 
in the case of domain wall motion dominated reversal is 
the field where domain wall motion starts depends on the 
size and on the shape of the nucleus. In an experimental 
situation point-like defects which are not visible in the 
optical regime can act as nuclei while there is no such 
artificial nucleus in our simulation. Second, some of the 
parameters of our model like the disorder or the Bloch- 
wall energy which influences the prefactors in Eqjl] may 
depend on the film thickness. In our simulation, they 
are thought to be constant - for simplicity and since we 
do not have more information about the thickness de- 
pendence of these parameters. Third, the time scale at 
which a hysteresis loop is observed may play an impor- 
tant role. This time scale is a few minutes in the exper- 
iments and around 60000 MCS in the simulations. To 
our knowledge, there is no straight forward way to map 
Monte Carlo simulation time on realistic time scales so 
that the observation-time windows might be different in 
experiment and simulation leading to different nucleation 
fields. Changing the sweeping rates changes the nucle- 
ation fields due to thermal activation, experimentally as 
well as in simulations, but this does not change the qual- 
itative behavior. 

However, it is not the aim of our simulations to calcu- 
late the nucleation field accurately. Rather, the simula- 
tions are thought to contribute to a better understanding 
of the fundamental properties of the system, especially to 
a better understanding of the dynamics. 

V. DYNAMICS AND TEMPERATURE 
DEPENDENCE 

The different reversal mechanisms we mentioned in the 
previous section manifest themselves also in a change of 
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the dynamical behavior. We focus on the reversal dynam- 
ics, i. e. the time dependence of the magnetisation after 
a rapid change of the field to a value which destabilizes 
the initial direction of magnetisation. The corresponding 
experimental and simulational time dependences of the 
magnetisation are compared in Fig. ^[ The time-axis is 
normalized to that time at which half of the system is 
reversed. 
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mains grows quadratic in time as long as the domains are 
so small that they do not overlap, M(t) — Mq ~ — (vt) 2 . 
Consequently, the central quantity of the magnetisation 
reversal driven by domain wall motion is the domain wall 
velocity. Therefore, in the following we will investigate in 
detail the dependence of the domain wall velocity on the 
driving field and the influence of the temperature. Hence 
from now on we restrict ourselves to the investigation of 
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FIG. 3. Domain image from MOKE measurements of a 
lOnm film 236, 313, and 468s after a rapid quench to a field 
of 17.7kA/m, a), and simulated domain configurations, of a 
150 x 150 system during the reversal of a lOnm film 10, 40, 
and 60 MCS after a rapid quench to a field of 245 kA/m, b) 



FIG. 2. Magnetisation versus time from measurement, a), 
after a rapid change of the field to 19.33kA/m (lOnm) and 
13.29kA/m (30nm) and from simulations, b), after a rapid 
change to 260 kA/m (lOnm) and 220 kA/m (30nm). 

For the case of the thin films the (absolute) value of 
the reversed field is larger than the coercive field while 
for the thicker film in the experimental situation a value 
above the nucleation field but below the coercive field 
had to be chosen since for larger fields the reversal is 
to fast to be observed with this experimental technique. 
Consequently, the long time limit of the corresponding 
experimental curve is above -1. Apart from that the ex- 
perimental results and the results from simulations agree. 

For the case of the 30nm film, i. e. for the nucleation 
driven reversal there is a rapid change of the magneti- 
sation at the beginning of the reversal process. This 
demonstrates that the reversal process is dominated by 
nucleation and not by domain growth processes. In the 
limit of a constant nucleation rate R and no growth pro- 
cesses the change of the magnetisation can be expected 
to be exponential p] . 

On the other hand for the case of domain growth the 
change of magnetisation is slow at the beginning. If a 
constant number of nuclei and a constant domain wall ve- 
locity v is assumed, the magnetisation of the reversed do- 



For the determination of the domain wall velocity 
within the simulation we start with a system that has 
a nucleus of circular shape with a radius of 19 cells in 
the center of the 150 x 150 system. When we switch 
on the driving field from the nucleus a domain starts to 
grow. For the better observation of the domain growth, 
in our flip-algorithm we do not consider cells that are 
not connected to the growing domain, i. e. we exclude 
the possibility of spontaneous nucleation. Otherwise we 
had - at least for finite temperatures - the problem that 
spontaneously new nuclei are build by thermal activa- 
tion which with increasing radius overlap with the origi- 
nal domain. Fig. |3| illustrates the time development of a 
domain which follows from the method described above 
and compares it with experimental domain images. The 
black regions are reversed domains following the mag- 
netic field. The domains look similar although their size 
is different (the linear size of the pictures is 3 /xm in the 
simulation and 150/^m experimental). This and also a 
more detailed analysis jl2| leads to the conclusion that 
the domain walls are rough and hence self-similar (see 
|i~3| for a review on self similar interfaces). The shape of 
the domains is - within the simulations - influenced by 
the strength of the disorder. Note, however, that frozen 
disorder is not the only possible reason for the rough- 
ening of a domain wall. Other possible sources are the 
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dipolar interactions as well as the dynamics of the model 
(see e.g. Jl3] for a discussion of the shape of domains in 
a model without disorder). 

From both, the domain configurations from simula- 
tions and from experimental domain images the mean 
radius r of the domains can be determined through the 
area F of the reversed domain as r = yj F/ir, assum- 
ing that the domain has a circular shape. Additionally, 
the domain radius of the experimental domain images 
has been determined as the mean distance of the domain 
boundary from the center of the domain. Both definitions 
of r lead to the same results. 

Fig. |^a shows the time dependence of the radius mea- 
sured at T — 299K after a rapid change of the field to 
H = 17.7, 18.0, and 18.8kA/m. Obviously, the domain 
wall velocity is constant and from the slope of the curves 
v can be determined. 
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FIG. 4. Measured radius of a domain versus time for 
H = 17.7,18.0,18.8 kA/m. T = 299K, a), and ra- 
dius of the domain versus time from simulations for 
H = 225, 230, 235, 240, 245 and zero temperature. The solid 
lines are best fitted. 

In order to get a deeper understanding of the influ- 
ence of the temperature on the dynamics apart from the 
"experimental temperatures" (T w 300K) we also simu- 
lated "extreme temperatures" (T = OK and T = 600K). 
Fig. [|b shows - as an example - the r(t) behavior from 
the simulations for T = OK and different fields. For the 
lowest field shown the domain wall is pinned, i. e. after 
a short period of rearrangement of the domain wall its 



movement stops and the radius remains constant. The 
pinning of the domain wall is due to energy barriers which 
follow from the disorder, the dipole field, and the intrinsic 
energy barrier of the single cell. For finite temperatures, 
the domain wall velocity is always finite (see discussion 
below). 

For r > 150L/2 = 1.5/xm the domain reaches the 
boundary of the system and - consequently - r(t) satu- 
rates. For smaller r the slope of the r(t) curve is approx- 
imately constant and v can be determined by fitting to a 
straight line. Fig. [5]shows the resulting dependence of the 
domain wall velocity on the driving field for T = 0, 300, 
and 600K. 
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FIG. 5. Simulated domain wall velocity versus driving field 
for T = 600,300, and OK (from above). The solid lines are 
guides to the eye. 

For zero temperature there is a sharp depinning tran- 
sition [ p"3[ at a critical field H c from a pinned phase with 
v = to a phase with finite domain wall velocity. This 
transition can be interpreted in terms of a dynamic phase 
transition with v ~ (H — H c ) 6 for H > H c where in our 
case the critical exponent is 9 ss 1, a value which is the 
mean field result for a moving elastic interface |l5f| in 
a random-field. Also, this value has been observed in 
simulations of a soft spin model with random-fields [^6| . 
Hence, it seems to be reasonable to consider the depin- 
ning transition we found to be in the universality class of 
a driven interface in the random-field Ising model. This 
fact is further confirmed by a dynamical scaling analysis 
of the structure of the domain walls in CoPt . The 
central quantity in this analysis is the time dependent 
roughness of the domain wall which - following the the- 
ory of driven interfaces Jl3[ - is characterized by a certain 
set of critical exponents. For the case of the CoPt alloy 
these exponents are also the same as those characterizing 
a driven interface in a random-field Ising model. 

Note that in the models mentioned above the only ori- 
gin for the depinning transition is the disorder. Here, the 
value of H c without disorder is zero. In our model, H c 
depends additionally on the dipole field and the intrinsic 
energy barrier of the single cell. A detailed analysis of 
these dependencies must be left for the future. 
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FIG. 6. Scaling plot from Fig. |. 

For finite temperatures the transition is smeared since 
for finite temperatures there is even for H < H c for each 
energy barrier a finite probability that the barrier can be 
overcome by thermal fluctuations. Hence, for finite tem- 
peratures we expect a crossover from the dynamics of the 
zero temperature depinning transition explained above to 
a dynamics dominated by thermal activation where the 
corresponding waiting times are exponentially large. For 
H < H c the domain wall velocity should decrease like 
Inv ~ (H - H c )/T. To illustrate this in Fig.| we show 
the corresponding semi-logarithmic scaling plot. As we 
expect, the data for the two different finite temperatures 
collapse for H < H c on a straight line. For H > H c ther- 
mal activation is obviously less relevant since for these 
values of the driving fields even for zero temperature the 
domain walls move. Consequently, also for finite temper- 
ature but H > H c the domain wall dynamics is domi- 
nated by the zero temperature depinning transition, i.e. 
v~H~H c . 
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FIG. 7. Scaling plot for the measured domain wall veloci- 
ties at T = 299, 335, and 394K. 

Fig. ^ shows the same scaling plot for the experimental 
data. Obviously, all data are in the regime H <C H c , so 
that the scaling for thermal activation works. We are far 
away from the regime H ps H c where the crossover to zero 
temperature depinning dynamics would start since from 
Fig. it follows H c ps 40kA/m. There is no possibility for 



measurements at higher fields since for higher fields new 
domains build up spontaneously and overlap each other 
so that it is not possible to follow one single domain for 
the analysis of the domain wall velocity. 

VI. CONCLUSIONS 

By comparing microscopy measurements on Co2sPt72 
alloy samples based on the magneto-optical Kerr effect 
with Monte Carlo simulations we show that a micromag- 
netic model can be used for the understanding of mag- 
netic properties of nanoscale magnetic films with high 
perpendicular anisotropy. The experimental and simula- 
tional results of the hysteresis, the reversal mechanism, 
the domain configurations during the reversal, the time 
dependence of the magnetisation and the temperature 
and field dependence of the domain wall velocity are in 
very good qualitative agreement. 

For thin films the reversal is dominated by a growth 
of domains the dynamics of which can be described by 
the domain wall velocity. The results for the domain 
wall velocity suggest that for zero temperature the hys- 
teresis can be understood as a depinning transition of 
the domain walls. For finite temperatures the transition 
is rounded by thermal activation and for fields smaller 
than the depinning field the domain wall movement is 
dominated by thermal activation. 
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